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tative signal of the global magnetic diurnal variation based on a spatial dis- 
tribution in both longitude and latitude of a set of magnetic stations as well 
as their magnetic behavior on a time basis. For that, we apply the Principal 
Component Analysis (PCA) technique implemented using gapped wavelet 
transform and wavelet correlation. The continuous gapped wavelet and the 
wavelet correlation techniques were used to describe the features of the mag- 
netic variations at Vassouras (Brazil) and other 12 magnetic stations spread 
around the terrestrial globe. The aim of this paper is to reconstruct the origi- 
nal geomagnetic data series of the H-component taking into account only the 
diurnal variations with periods of 24 hours on geomagnetically quiet days. 
With the developed work, we advance a proposal to reconstruct the baseline 
for the quiet day variations (Sq) from the PCA using the correlation wavelet 
method to determine the global variation of PCA first mode. The results 
showed that this goal was reached and encourage other uses of this approach 
to different kinds of analysis. 

Keywords: Magnetogram data, H-component, Quiet days, Principal 
Component Wavelet analysis. 



i 1. Introduction 



2 A substantial part of the energy carried by the solar wind can be trans- 

3 fered into the terrestrial magnetosphere and it is associated with the passage 



of southward directed interplane tary magnetic fi e 



5 ficiently long intervals of time. 



Gonzalez et ai- 



ds, B s, by the Earth for suf- 



( 119941 ) discussed the energy 



e transfer process as a conversion of the directed mechanical energy from the 
7 solar wind into magnetic energy stored in the magnetotail of Earth's magne- 
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8 tosphere and its reconversion into thermal mechanical energy in the plasma 

9 sheet, auroral particles, ring current, and Joule heating of the ionosphere. 

10 The increase on the solar wind pressure is responsible for the energy in- 

11 jections and induces global effects in the magnetosphere called geomagnetic 

12 storms. The characteristic signature of geomagnetic storms can be described 

13 as a depression on the horizontal component of the Earth's magnetic field 

14 measured at low and middle latitude ground stations. The decrease in the 

15 magnetic horizontal field component is due to an enhancement of the trapped 

16 magnetospheric particle population, consequently an enhanced ring of cur- 

17 rent. This perturbation of the H- component could last from several hours to 



is several days (as described by 



Kamide et al. 



1998|). 



19 The geomagnetic storms can consist of four phases: sudden commence- 

20 ment, initial phase, main phase and recovery phase. The sudden commence- 

21 ment when it exists, corresponds to the moment when the initial impact of 

22 the increased solar wind pressure over the magnetopause occurs. The ini- 

23 tial phase at ground appears as a rapid increase on the H-component over 

24 less than 1 h almost simultaneously worldwide. The main phase of the ge- 

25 omagnetic storm lasts a few hours and is characterized by a decrease in the 

26 H-component. The recovery time corresponds to the gradual increase of the 



27 H-component value to its average level. A de tailed description of t 



ne mor- 



28 phology of magnetic storms is, for instance, in I Gonzalez et al.l (jl994j ). 

29 The intensity of the geomagnetic disturbance in each day is described 

30 by indices. The indices are very useful to provide the global diagnostic of 

31 the degree of disturbance level. There are different indices that can be used 

32 depending on the character and the latitude influences in focus. Considering 
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only the main latitudinal contributions, the ring current dominates at low 
and middle latitud es and an auroral iono spheric current systems dominates 



35 at higher latitudes (IMendes et al 



2005a] ). Kp, AE and Dst and their deriva- 



tions are the most used geomagnetic indices. The Kp index is obtained from 
the H-component and it is divided in ten levels from to 9 corresponding to 
the mean value of the disturbance levels within 3-h intervals observed at 13 
subauroral magnetic stations (see 



B artels 



19571 ). However, the K index is 



the most difficult to be physically interpreted due to its variations be caused 
by any geophysical current system including magnetopause currents , field- 



aligned currents, and the auroral electrojets (I Gonzalez et al. 



199J). The 



minutely AE index (sometimes 2.5 minute interval) is also obtained by the 
H-component measured from magnetic stations (5 to 11 in number) located 
at auroral zones and widely distributed in longitude. The AE index provides 
a measure of the overall horizontal auroral oval current strength. The index 
most used in low and mid-latitudes is the Dst index. It represents the varia- 
tions of the H-component due to changes of the ring current and is calculated 
every hour. 

The Dst index is described as a measure of the worldwide derivation 
of the H-component at mid-latitude ground stations from their quiet days 
values. At mid-latitude, the H-component is a function of t he ma gnetopause 



Burton et al. 



(119751 ) calculated 



currents, the ring current and tail currents, 
the Dst index as a average of the records from N mid-latitude magnetic 
stations following, 



Dst 



1 N 

i=l 



disturbed 



-H, 



quiet 



AH 



(1) 
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where AH is a local time H average, Hdisturbed is the H-component mea- 
sured at disturbed days and H quiet , on quiet days. 

Other contributions beyond th e ring current c ould be extracted or elim- 



59 inated with the idea presented by 



Burton et al. 



(119751 ). Those authors de- 



60 scribed the evolution of the ring current by a simple first order differential 

ei equation, 

d Dst* 

a ^ = Q{t)-aDst\ (2) 

62 where Dst* = Dst — bPdyn^ + c. The contribution of the magnetopause 

63 currents to H is proportional to the square root of the solar wind dynamic 

64 pressure (Pdyn), Q represents the injection of particles to the ring current, 

65 a Dst* represents the loss of particles with an e-folding time - and the con- 
ee stant terms a, b and c are determine by the quiet days values of the magne- 
67 topause and ring currents. 

6B The Dst index is available on the Kyoto World Data Center at http:// 

69 wdc.kugi.kyoto-u.ac.jp/dstdir/index.html. It is traditionally calculated from 

70 four magnetic observatories: Hermanus, Kakioka, Honolulu, and San Juan. 

71 These observatories are located at latitudes below 40° which are sufficiently 

72 distant from the auroral electrojets. The derivation of the Dst index corre- 

73 sponds to three main steps: the removal of the secular variation, the elimi- 

74 nation of the Sq variation and the calculation of the hourly equatorial Dst 



75 Index (see http://wdc.kugi.kyoto-u.ac.jp/dstdir/dst2/onDstindex.html) 



76 The traditional method of calculating the baseline for the quiet day vari- 

77 ations uses the five quietest day for each month for each magnetic observa- 
7B tory. In this work, we propose a way to deal with Sq variations by suggesting 
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a method using Principal Component with the wavelet correlation matrix. 
This method eliminates the disturbed days using a multiscale process. Also, 
we developed an algorithm for extracting the solar quiet variations recorded 
in the magnetic stations time series, in order words, a way of estimation of 
the quiet-time baseline. To accomplish this task, we separate the solar di- 
urnal variatio ns using hourly data of the H-component using the technique 



85 (described in 



Klausner et al. 



2011aj). Afterward we applied the principal 



component wavelet analysis to identify the global patterns in the solar diur- 
nal variations. 

The rest of the paper is organized as follows: Section [2] is devoted to 
explain the main issues of the Dst index calculation procedure. In Section [3] 
the analyzed period and data are presented. Section H] describes the principal 
component analysis and it is devoted to introduce the suggested method of 
principal component analysis (PCA) using gapped wavelet transform and 
wavelet correlation. It also establishes the identification of the disturbed 
days. The results are discussed in Section [51 and finally, Section [6] brings the 
conclusions of this work. 



2. The Dst index calculation procedure 



Recently, reconstructing the Dst a nd removing the quiet-time baseline 



98 have been a motivation of several wo rks (K arinen and Mursulal. 



2005a; 



2009 



Karinen and Mursula 



Klausner et al 



2006 



Mursula et al 



2008 



2005 



Mendes et al. 



Love and Gannon . 



2011al ). Despite of these disagreements in the Dst con- 

101 strution, today, the Dst index remains an important tool in the space weather 

102 analysis. 
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Karinen and Mursulal ( 120051 ) reconstructed the Dst index (Dxt) follow- 
ing the original formula presented at http:/ /wdc.kugi.kyoto-u.ac.jp| /dst- 
dir/dst2/onDstindex.html. However, they encountered a few issues as: the 
availability and the data quality , some shifts in the baseline level of the 
H-component, erroneous data points and some data gaps at all/some mag- 
netic observatories. The Dst could not be fully reproduced using the original 
formula because the inadequate information above effects the treatment of 
the related issues and there fore remains partly unscientific as described by 



Karinen and Mursulal (120051 ). 



A ne w corrected and extended ve rsion (Dcx) of the Dxt index was pro- 



113 posed by 



Karinen and Mursulal (120061 ) . They corrected the Dst index for the 



excessive seasonal varying quiet-ti me level which was un r elated to magnetic 
storms as previously discussed in iKarinen and Mursulal (120051 ) . They also 
showed that the considerable amount of quiet-time variation is included in 
the Dxt index but none in the Dcx index. 

Another issue related to the derivation of the Dst index is that no treat- 



ment is made 
observatories. 



:o normalize t 



Mursula et ai- 



r e diffe rent latitudinal location of the magnetic 
(120081 ) suggested the normalization of the mag- 
netic disturbances at the four Dst stations with different latitudes by the co- 
sine of the geomagnetic latitude of the respective station. If no correction is 
made, they showed for the lowest geomagnetic station, Honolulu, the largest 
deviations and for the highest station, Hermanus, the lowest deviations of 
the four station. The standard deviations reflect the annually averaged ef- 
fect of the (mainly ring current related) disturbances at each station, (see 



Mursula et al. 



20081 . for more details). 
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Mendes et al.l (j2005al ) evaluated the effect of using more than four mag- 
netic stations and shorter time intervals to calculate the Dst index. The 
obtained Dst index profiles using 12, 6 or 4 magnetic station did not show 
significant discrepancies and the best agreement with the standart Dst was 
obtained using magnetic stations located at latitudes lower than 35° in both 
hemispheres. 

Although, the increase of symmetrically world-wide distributed magnetic 
stations did not effect significantly the Dst index, the longitudinal asymme- 
tries of the ring current contributes for t he average disturbanc es of the Four 



Mursula et al. 



(120101 ) using an ex- 



Dst stations be systematically different, 
tended network of 17 stations, demonstrated that the local disturbances are 
ordered according to the station's geographic longitude, where the western- 
most station (Honolulu) presented the largest disturbances and contributions 
to Dst index and the east ernmost (Kakioka) the smallest. 

Klausner et al.l (j2011al ) studied the characteristics of the Sq variations at 
a Brazilian station and compared to the features from other magnetic sta- 
tions to better understand the dynamics of the diurnal variations involved in 
the monitoring of the Earth's magnetic field. They used gapped wavelet anal- 
ysis and the wavelet cross-correlation technique to verify the latitudinal and 
lo ngitudinal dependence of the diurnal variation s. As previously mentioned 



by 



Mursula et al. 



(120101 ) 



Klausner et al. 



( 1201 laf ) also verified that magnetic 



stations located at lower latitudes and westernmost (Honolulu and San Juan) 
presented larger correlation to Vassouras than the easternmost stations (as 
Kakioka) . 

Some important aspects for the construction of the Dst index as described 
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153 by lLove and Gannonl (120091 ) are: the utilization of the original data, the in- 

154 spection in time and frequency domains (removal of diurnal variation) and 

155 the consideration of the distinction between stationary and non-stationary 

156 time series ingredients which applie s to the secular variation. Also as men- 



157 tioned by 



Love and Gannon 



(120091 ). some patterns of the global magnetic 

158 disturbance field are well understood and some are not which means that 

159 there is still a lot to learn about the magnetosphere, magnetic storm and 
wo Earth-Sun relationship. 



161 3. Magnetic Data 

162 In this paper, we use ground magnetic measurements to estimate the 

163 quiet-time baseline. We select the four magnetic observatories used to cal- 

164 culate the Dst index: Hermanus (HER), Kakioka (KAK), Honolulu (HON), 

165 and San Juan (SJG), plus other 9 different magnetic observatories reasonably 
lee homogeneously distributed world wide. One of these nine chosen stations is 
167 Vassouras (VSS) located under the South Atlantic Magnetic Anomaly (min- 
us imum of the geomagnetic field intensity). The geomagnetic data use in this 
leg work relied on data collections provided by the INTERMAGNET programme 
no (htt p:/ /www.intermagnet.org[ ). 

ni The distribution of the magnetic stations, with their IAGA codes, is given 

172 in Fig. [TJ The corresponding codes and locations are given in Table [TJ These 

173 selection of mag netic stations correspo nd to the same selection used in a 

174 previous work o felausner et al.l . l2011al ) for the same reasons (exclusion of 

175 the major influence of the auroral and equatorial electrojets). In this work, 
ne we only use the data interval corresponding to the year 2007. We also ap- 
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177 ply the same methodo logy to identify geomagnetically quiet days used by 



Klausner et al. 



(1201 lal ) . We consider quiet days, only those days in which 
179 the Kp index is not higher than 3+. 

wo As at low latitudes the horizontal component (H) is mostly affected by 
lsi the intensity of the ring current, we decided to use only the hourly mean 

182 value series of this component. The magnetic stations present available data 

183 in Cartesian components (XYZ system). The Co nversion to horizo ntal-polar 



1997, for more 



184 components (HDZ system) is very simple (see ICampbelll . 

185 details). The system's conversion was performed in all the chosen magnetic 

186 stations. 




-60 



-150 -120 



120 150 



Figure 1: Geographical localization of the stations used in this work and their respective 
IAGA code. 



is? 4. Methodology 

188 The method used in this study is based on the principal component analy- 

189 sis (PCA) using gapped wavelet transform and wavelet correlation to charac- 
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Table 1: INTERMAGNET network of geomagnetic stations used in this study. 



Station 


Geographic coord. 


Geomagnetic coord. 




Lat.(°) 


Long.(°) 


Lat.(°) 


Long. ("J 


AMS 


-37.83 


77.56 


-46.07 


144.94 


ASP 


-23.76 


133.88 


-32.50 


-151.45 


BEL 


51.83 


20.80 


50.05 


105.18 


BMT 


40.30 


116.20 


30.22 


-172.55 


BOU 


40.13 


-105.23 


48.05 


-38.67 


CLF 


48.02 


2.26 


49.56 


85.72 


CMO 


64.87 


-147.86 


65.36 


-97.23 


EYR 


-43.42 


172.35 


-46.79 


-106.06 


HER 


-34.41 


19.23 


-33.89 


84.68 


HON 


21.32 


-158.00 


21.59 


-89.70 


KAK 


36.23 


140.18 


27.46 


-150.78 


SJG 


18.12 


-66.15 


27.93 


6.53 


vss 


-22.40 


-43.65 


-13.43 


27.06 



Source: http:/ /wdc. kugi.kyoto-u.ac.jp/igrf/gggm/index. html (2010) 



wo terize the global diurnal variation behavior. To identify periods of magnetic 

191 disturbance, we use the discrete wavelet transform. This technique is em- 

192 ployed to analyze the removal of disturbed days from the magnetograms, and 

193 consequently, from the reconstructed Sq signal. Also in this section, a com- 
w bined methodology using the PCA and gapped wavelet transform is briefly 
195 described. Following, we present an identification method to distinguish the 
we disturbed days using discrete wavelet coefficients. 
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Global geomagnetic behavior analysis 

Among the several available methods of analysis, PCA is a particularly 
useful tool in studying large quantities of multi-variate data. PCA is used 
to decompose a time-series into its orthogonal component modes, the first of 
which can be u s ed to describe the dominant patterns of variance in the time 



202 series 



Murray 



( 1l984j ). The PCA is able also to reduce the original data set 
of two or more observed variables by identifying the significant information 
from the data. Principal Components (PCs) are derived as the eigenvectors of 
the correlation matrix between the variables. Their forms depend directly on 
the interrelationships existing within the data itself. The first PC is a linear 
combination of the original variables, which when used as a linear predictor 
of these variables, explains the largest fraction of the total variance. The 
second, third P C, etc., explain the largest parts of the remaining variance 



Murray 



fll984h . 



211 As explained by lYamadal (l2002f ). the interpretation of the eigenvectors 

212 and the eigenvalues can be described as follow, the eigenvectors are the nor- 

213 malized orthogonal basis in phase space, and also, the set of vectors of the 

214 new coordinate system in the space, different from the coordinate system of 

215 the original variables; the eigenvalues are the corresponding variance of the 

216 distribution of the projections in the new basis. 

217 In order to isolate the global contributions of each PCs mode, we applied 

218 PCA using the wavelet correlation matrix computed by ga pped wavelet trans- 



form. This wavelet correlation matrix was introduced in 



Nesme-Ribes et al 



220 ( 119951 ) . We joined the properties of the PCA, which are the compression of 



large databases and the simplification by the PCs modes, and properties of 
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222 the wavelet correlation matrix, which is the correlation at a given scale, a, 

223 in this case, the scale corresponded to the pseudo-period of 24 hours. 

224 Identification of magnetic disturbance 

225 The wavelet analysis has the following propriety: the larger amplitudes of 

226 the wavelet coefficients are associated with l ocally abrupt s i gnal ch anges or 



"details" of highe r frequency. In the wor k of 



228 following work of iMendes da Costa et al.l (120111 ). a method for the detection 



Mendes et al. 



(J2005bJ) and the 



229 of the transition region and the exactly location of this discontinuities due 

230 to geomagnetic storms was implemented. In these cases, the highest ampli- 

231 tudes of the wavelet coefficients indicate the singularities on the geomagnetic 

232 signal in association with the disturbed periods. On the other hand, when 

233 the magnetosphere is under quiet conditions for the geomagnetic signal, the 

234 wavelet coefficients show very small amplitudes. In this work, we applied this 

235 methodology with Daubechies orthogonal wavelet function of order 2 on the 

236 one minute time resolution with the pseudo-periods of the first three levels 

237 of 3, 6 and 12 minutes. 

238 5. Results and Discussion 

239 In this section, we will present the results of reconstructed baseline for the 

240 global quiet days variation using PCA technique implemented with gapped 

241 wavelet transform and wavelet correlation. Also, we will apply DWT to eval- 

242 uate the day-by-day level of geomagnetic disturbance using KAK magnetic 

243 station as reference. 

244 Fig. [2] shows an example of the geomagnetic behavior presented at the 

245 June, 2007 magnetogram of Kakioka using continuous gapped wavelet trans- 
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Figure 2: The GWT analysis of KAK magnetogram. At top, it shows the H-component 
of KAK at June, 2007 used for the wavelet analysis, at bottom left, the scalogram using 
Morlet wavelet, logarithmic scaled representing log2(|T / V r (o,6)|), and at bottom right, the 
global wavelet spectrum. 

246 form (GWT). The GWT can be used in the analysis of non-stationary signal 

247 to obtain information on the frequency or scale vari ations and to detect its 



Klausner et al 



2011a 



248 structures localization in time and/or in space (see 

249 for more details). It is possible to analyze a signal in a time-scale plane, 

250 called so the wavelet scalogram. In analogy with the Fourier analysis, the 

251 square modulus of the wavelet coefficient, |W(a, b)\ 2 , is used to provide the 

252 energy distribution in the time-scale plane. In the GWT analysis, we can also 

253 explore the central frequencies of the time series through the global wavelet 

254 spectrum which is the variance average at each scale over the whole time 

255 series, to compare the spectral power at different scales. This figure shows 

256 the H-component (top), the wavelet square modulus (bottom left) and the 

257 global wavelet spectrum (total energy in each scale - bottom right). In the 
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258 scalogram, areas of stronger wavelet power are shown in dark red on a plot 

259 of time (horizontally) and time scale (vertically). The areas of low wavelet 

260 power are shown in dark blue. 

261 In Fig. |2l it is possible to notice peaks of wavelet power on the scalogram 

262 at the time scale corresponding to 8 to 16 minutes of period. This periods are 

263 associated to PC5 pulsations during disturbed periods. Also, it is possible 

264 to notice a maximum of wavelet power at the time scale corresponding to 

265 harmonic periods of the 24 hours such as 6, 8, 12 hours. Those periods are 

266 related to the diurnal variations. 

267 The GWT technique is able to analyze all the informations present on the 

268 magnetogr ams. 



t is an auxiliary tool to localize on time/space the PC1-PC5 

269 pulsations ( ISaitol . Il969h . However, the scalogram provides a very redundancy 

270 information which difficult the analysis of each decoupling phenomenum. For 

271 that reason, we preferred to use DWT to evaluate the day-by-day level of 

272 geomagnetic disturbance. 

273 Fig. [3] is composed of two graphs, the first one presents the reconstructed 

274 line of the Sq variation where the first 10 geomagnetically quietest days of 

275 each month are highlighted in blue and the 5 most disturbed days in red 

276 and the second one presents the discrete wavelet analysis of the geomagnetic 

277 horizontal component obtained at Kakioka station, Japan. 

278 Using our criteria of removing disturbed days, we consider as gaps of 3rd, 

279 14th, 21th and 29th day. In our case, the gapped wavelet technique is very 

280 helpful because it reduces two effects: the presence of gaps an d the boundary 



281 effects due to the finite length of the data, for more details see 

282 faplli . 



Klausner et al 
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Figure 3: The amplitude variation of the reconstructed Sq signal with the highlighted 10 
geomagnetically quietest days (blue) and 5 most disturbed days (red), the H-component 
average variation for KAK obtained at June, 2007 and the square root wavelet coefficients 
amplitudes dl, d2 and d3 with the pseudo-periods of 3, 6 and 12 minutes. 

283 The first graph shows the amplitude range between —5 and 5 nT and 

284 presents a complex pattern. It is possible to notice that the larger amplitudes 

285 of the reconstructed Sq signal correspond to the periods between the days 

286 8-10, 13-15, 21-23 and 29-30. These periods correspond to the disturbed 

287 days. 

288 Table [5] shows these 10 quietest days and 5 most disturbed days of each 

289 month set by the GeoForschungsZentrum ( GFZ) Potsdam through the anal- 

290 ysis of the Kp index that are highlighted on the second graph. The year 

291 of 2007 is a representative year of minimum solar activity and it is used in 

292 this analysis due to have less disturbed periods (see our considerations in 

293 Section [3]). By analyzing these highlighted days, we expect to find out if 

294 there is a correlation between the days classified as quiet days and a small 
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295 Sq amplitude variation. 



Table 2: The first 10 geomagnetically quietest days and first 5 most disturbed days set by 
the GeoForschungsZentrum ( GFZ) Potsdam 



Month/year 






10 quietest days 








5 most disturbed days 




ql 


q2 


q3 


q4 


q5 


q6 


q7 


q8 


q9 


qlO 


dl 


d2 


d3 


d4 


d5 


Mar 


20 


21 


3 


19 


9 


29 


22 


18 


31 


4 


13 


24 


6 


7 


14 


Jun 


5 


12 


6 


7 


11 


26 


25 


20 


30 


28 


14 


21 


22 


3 


29 


Sep 


13 


9 


10 


11 


12 


16 


17 


19 


26 


15 


29 


2 


28 


23 


27 


Dec 


3 


8 


4 


25 


7 


26 


2 


29 


15 


6 


18 


17 


11 


20 


21 



SOURCE: http://wdc.kugi.kyoto-u.ac.jp/qddays/index.html 

The second graph shows the discrete wavelet analysis applied to geo- 
magnetic minutely signal from KAK using Daubechies orthogonal wavelet 
family 2. From top to bottom in this graph, the H-component of the ge- 
omagnetic field and the first three levels of the square wavelet coefficients 
de noted by dl, d2 and d3. This analysis uses the methodology developed 



Mendes et al. 



302 (1201 If ) and 



(200 5M), and pos 



Klausner et al. 



teriorly applied by 



Mendes da Costa et al 



(12011bD . 

In order to facilitate the evaluation of the quiet periods obtained by the 
discrete wavelet analysis applied to geomagnetic signal from KAK, we also 
developed a methodology (effectiveness wavelet coefficients (EWC)) to in- 
terpret the results shown in Fig. HJ The EWC corresponds to the weighted 
geometric mean of the square wavelet coefficients per hour. It is accom- 
plished by weighting the square wavelet coefficients means in each level of 



309 decomposition as following 

4 Ef=i dl + 2 Ef=i d2 + Ef =1 d3 



EWC 



7 



(3) 



310 where iV is equal to 60 because our time series has one minute resolution. 
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Figure 4: The comparative of the global Sq behavior and the effectiveness wavelet coeffi- 
cients for the month of June, 2007. 

311 Through Fig. [U it is possible to compare the global Sq behavior (top 

312 graph) with the analysis of quiet and disturbed days obtained by one repre- 

313 sentative magnetic station of medium/low latitudes (KAK - bottom graph) 
3w in order to verify situations in which the global Sq behavior presents less 

315 or more variability. This analysis allows us to validate the quietest days, 

316 and evaluate the most disturbed days in order to establish a reliable method 

317 of global Sq analysis obtained from medium/low latitude magnetic stations 

318 influenced only by the ionosphere. 

319 In Fig. HI the global Sq behavior (top graph) shows larger amplitudes 

320 during the periods between the days 8-10, 12-16, 22-23 and 29-30. Most 

321 of these periods correspond to the days where the EWCs have an increase 

322 of their values as shown KAK analysis (bottom graph). The increase of the 

323 EWCs values occurs during the periods between the days 1-4, 8-10, 13-18, 
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324 21-24 and 29-30. The EWCs can help us also to interpret the results obtained 

325 in each day, and can help us to evaluate the quietest and most disturbed days 

326 measure by the selected magnetic station of medium/low latitudes, KAK, as 

327 shown in Table |5j 

Table 3: The 10 geomagnetically quietest days and 5 most disturbed days obtained by the 



discrete wavelet analysis 



Month/year 
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5 most disturbed days 
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q5 


q6 
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13 
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12 
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Jun 


11 
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25 
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26 
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19 
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21 
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14 


22 


29 
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12 


11 
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17 
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16 


13 


26 


19 


27 


22 


29 
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20 
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25 
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1 


21 
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6 


16 


28 


17 


18 


11 


20 


10 



328 The same methodology and analysis comparing the global Sq behavior 

329 and the EWCs from KAK is done for the month of March, September and 

330 December, 2007. 

331 Fig. shows the comparative of the global Sq behavior and the EWCs 

332 from KAK for the month of March, 2007. The amplitude range of global 

333 Sq signal is between —15 and 15 nT, and, it also shows a complex pattern. 

334 It is possible to notice that the larger amplitudes of the reconstructed Sq 

335 signal correspond to the periods between the days 6-7, 11-17 and 23-28. 

336 Once more, these periods correspond to the disturbed days. The increase of 

337 the EWCs values occurs during the periods between the days 1, 6-7, 11-17, 

338 23-28 and 30-31. We can notice that the increase of the amplitude of the 

339 reconstructed Sq signal correspond to the increase of the EWCs magnitude. 

340 The amplitude range for September, 2007, is between —15 and 15 nT 

341 and the larger amplitudes of the reconstructed Sq signal correspond to the 

342 periods between the days 1-7, 19-23 and 26-30, see Fig, [6J The increase of 

19 



Time {days) 



Figure 5: The comparative of the global Sq behavior and the effectiveness wavelet coeffi- 
cients for the month of March, 2007 

343 the EWCs values occurs during the periods between the days 1-3, 20-24 and 

344 27-30. Comparing these two analysis, the global and the EWCs, we verify 

345 that the KAK magnetic behavior represents well the increase of global Sq 

346 oscillations. 

347 In Fig. [7J the amplitude range is between —10 and 10 nT and the larger 

348 amplitudes correspond to the periods between the days 15-23. Also, the 

349 increase of the EWCs values occur during the periods between the days 10- 

350 11 and 17-21. Once more, when we compare these two analysis, the global 

351 and the EWCs, we verify that the KAK magnetic behavior represents well 

352 the increase of global Sq oscillations. 

353 We observe in Figs. HI [5j [6] and UJ in most of the cases, the major am- 

354 plitude fluctuations of the reconstructed Sq signal correspond to the most 

355 disturbed days and minor fluctuations, to the quietest days. However, the 
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Figure 6: The comparative of the global Sq behavior and the effectiveness wavelet coeffi- 
cients for the month of September, 2007 

356 diurnal global variability shows a complexity on the amplitude variation pat- 

357 tern even during geomagnetically quiet periods. Through this study we com- 

358 pare the amplitude variation of the reconstructed Sq signal to effectiveness 

359 wavelet coefficients obtained at KAK with the purpose of understanding the 

360 complexity of the diurnal global variability 

361 6. Conclusions 

362 In this work, we suggest an alternative approach for the calculation of the 

363 Sq baseline using wavelet and PCA techniques. This new approach address 

364 some issues, such as, the availability and the quality of data, abrupt changes 

365 in the level of the H-component, erroneous points in the database and the 

366 presence of gaps in almost all the magnetic observatories. To fulfill this 

367 purpose, we reconstruct the Sq baseline using the wavelet correlation matrix 
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Figure 7: The comparative of the global Sq behavior and the effectiveness wavelet coeffi- 
cients for the month of December, 2007 

368 with scale of 24 hours (pseudo-period). The PCA/wavelet method uses the 

369 global variation of first PCA mode that also corresponds to phenomena with 

370 periods of 24 hours. This study shows that the largest amplitude oscillation 

371 of the reconstructed signal (Sq baseline) corresponded to the most disturbed 

372 days and the smaller oscillations to the quietest days. This result is consistent 

373 with the expected Sq variations. 
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